Exactly solvable model of electron in the Lame potential 
and singularities of the electron thermodynamic potential 
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One-gap and two-gap separable Lame potentials are studied in detail. For the one-dimensional 
case, we construct the dispersion relation graph E{k) and for the three-dimensional case we construct 
the Fermi surfaces in the first and second bands. The pictures illustrate a passage from the limit 
case of free electrons to the limit case of tight binding electrons. These results are used to describe 
the Lifshits electron phase transition of 2| kind and derive some exact expressions. We also examine 
the singularities of the second derivative of magnetic momentum in an external magnetic field. The 
parameter of the singularities depends on corresponding effective mass. 
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I. INTRODUCTION 

There is a sufficient number of numerical and self- 
consistent methods to compute the one particle spectrum 
of a solid with good accuracy. However, in addition to 
this, it is desirable to have analytical solutions of the 
Schrodinger equation. These solutions for multidimen- 
sional periodic potentials are useful both for a zero order 
perturbation theory and for the general theory. This al- 
lows us, for example, to gain a deeper understanding of 
spectral manifolds and corresponding eigenfunctions, as 
well as their geometrical and analytical properties. 

Among onc-dimensional potentials, the finite-gap ones 
are the most interesting and fruitful in the view of the 
physical results. The Schrodinger equation with the 
finite-gap potential^i^i has an exact solution, the spec- 
trum and eigenfunctions are defined by analytical func- 
tions. 

The Kronig-Penny potentiaP^ consisting of rectangu- 
lar wells is usually considered as the most elementary 
one in the quantum physics of solids. However, in order 
to determine the band structure that corresponds to this 
potential, it is necessary to solve transcendental equation 
and, furthermore, the appropriate eigenfunctions are too 
cumbersome to use them to calculate the matrix elements 
of any observables. The finite-gap potentials are free of 
these limitations. The matrix elements of any observ- 
able on the eigenfunctions that correspond to finite-gap 
potential can be calculated analytically by the residue 
method. Furthermore, one can say that these potentials 
play the same role in solid-state physics as the Kepler 
problem plays in atomic theory. 

It is known that generic periodic potential has gen- 
erally an infinite set of gaps in the spectrum, but their 
width rapidly decreases as energy increases. If a poten- 
tial is analytical, then the rate of decrease is exponential. 
Therefore, provided narrow enough gaps in the spectrum 
are disregarded, any periodic potential can be approxi- 
mated by a finite-gap potential.-- It turns out also that 
the finite-gap potentials are exact solutions of the Pcicrls- 
Frohlich problem on a self-consistent state of the lattice 
and conduction electronsi^ 



The simplest finite-gap potentials are the Lame 
They are expressed through the Weier- 
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strass function. The Lame potentials are a subset of the 
Darboux potentials^ that are linear superpositions of 
the Weierstrass functions with half-period shifts. What- 
ever the number of gaps, all these potentials are char- 
acterized by only two parameters of the elliptic func- 
tion. Treibich and Verdier rediscovered the Darboux 
potentials and described them in algebraic-geometric 
terms of the tangent covers and studied their spectral 
properties Krichever suggested considering ellip- 

tic finite-gap potentials from the point of view a torus 
covering.— Using this approach, Smirnov has presented 
a large list of these potentials. Its and Matveev de- 
rived a general formula that expresses an arbitrary finite- 
gap potential through the multidimensional Riemann 9- 
functions, the parameter number of which coincides with 
the number of band edges.— 

There are also multidimensional elliptic Calogero- 
Moser potentials, which their integrability was proved 
by Olshanetsky and Perclomovi^ 

Baryakhtar, Belokolos and Korostil suggested consid- 
ering the Schrodinger equation with three-dimensional 
separable potential as a sum of three one-dimensional 
finite-gap potentials along orthogonal directions.^ It is 
obvious that these potentials correspond to the or- 
thorhombic lattice. The Schrodinger equation with such 
a potential has the exact eigenfunctions and the spectrum 
is defined by analytical functions. The parameters that 
this model includes have simple physical sense. These 
parameters are related to the width of spectral gaps. As 
an example, let us consider separable elliptic finite-gap 
potentials. The Weierstrass function has two indepen- 
dent parameters: the parameter uj' defines the poten- 
tial period, the parameter lu defines the gap width along 
corresponding direction of dual lattice. It is natural to 
consider the parameter t = lo' /lo as a characteristic of 
the relative width of gaps. As the value of |t| increases, 
the corresponding gap width increases too. In one of the 
extreme cases, the width of all the gaps vanishes. It cor- 
responds to the free electron model, while in the other 
extreme case, the width of all the bands vanishes and this 
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corresponds to the tight binding model. In the model of 
separable finite-gap potentials, we can also easily study 
one- and two-dimensional lattices. This can be done by 
choosing extremely wide gaps along certain directions. 

Finite-gap potentials have successfully been applied to 
solve different problems in solid-state physics.— '^'^lii'^ii^ 

In our paper, on one hand, we examine the known one- 
gap and two-gap separable Lame potentials in detail. For 
the one-dimensional case, we construct the dispersion re- 
lation graph E{k) and for the three-dimensional case, 
we construct the Fermi surfaces in the first and second 
bands. The following pictures illustrate a change from 
one extreme case of free electrons to the other extreme 
case of tight binding electrons. On the other hand, we 
use the separable one-gap Lame potential to describe the 
Lifshits electron phase transition of 2^ kind and derive 
some exact expressions. At the end, we consider singular- 
ities of the second derivative of the magnetic momentum 
of a normal metal in an external magnetic field. 

The paper is organized as follows. In section 2, we 
present information on the one-dimensional Lame poten- 
tials using two different approaches and the one-gap and 
two-gap ones are studied in detail. In section 3, we con- 
sider the separable finite-gap Lame potentials and con- 
struct the Fermi surfaces using the one-gap one. In sec- 
tion 4, the one-gap separable Lame potential is applied 
to describe the Lifshits electron phase transition of 2^ 
kind. At the end of this section we consider singularities 
of the second derivative of the magnetic momentum of a 
normal metal in an external magnetic field. 

II. THE ONE-DIMENSIONAL FINITE-GAP 
LAME POTENTIALS 

The Lame potentials can be represented in one of the 
two forms (since these potentials are expressed through 
the elliptic functions which have two independent peri- 
ods) 

vix) = eie + i)p{x + tj'\uj,u}'), (1) 

U{x) ^-i{e+l)p{ix + Lj\LJ,Lj'), (2) 

were ^ is a number of gaps in the spectrum, and 
p {z\uj, uj' ) is the elliptic Weierstrass function with a real 
half-period u) and an imaginary one lu' . We use the stan- 
dard notation of the elliptic function theoryi^ If it doesn't 
cause misunderstanding we will write simply p (z). 

The potential V{x) has the period 2uj, the period of 
the potential U{x) is 2\llj'\. Since the Weierstrass func- 
tion is homogeneous and it isn't changed under uni- 
modal transformations, the potentials H]) and ^ are ex- 
pressed one through the other: i{i + l)p{x + Cj') — 
-i{£ + l)p{ix + iu;'\ioj, iuj' ) = -£{£ + l)p{ix + uj' \uj, uj'). 
We put UJ = \uj'\, and uj' — iCo; correspondingly 77 = 
-\f)'\, 77' = -if] and 53 = -53, = 92, also d = 
—63, 63 = — ei, 62 = —62. In order to pass from the 



spectrum of the Schrodinger operator with the potential 
U{x) to the spectrum of the Schrodinger operator with 
the potential V{x), it is necessary to pass from uj, uj' to 
UJ, uj' in all the expressions. This means that we have 
to write —E instead of E, rj' /uj' instead of rj/uj and ik 
instead of k in expressions for the spectrum that will be 
defined later. 

We will examine the potentials U{x) only. 

The dimensionaless Schrodinger equatioit^l 

- d^-^ - £ {£ + 1) p {ix + uj) = E'i' (3) 

has two Bloch eigenfunctions which correspond to differ- 
ent signs of the wavevector 

In order to pass to dimensional values it is necessary to 
multiply E, U{x) by Eq — T? j^ma} and spacial variable 
X by a, were a is a characteristic length of the crystal 
potential; m is a particle mass and h is the Planck con- 
stant. 

W(E) is the Wronskian of the functions '^-(x,E), 
^j^(x,E), and W'^(E) is a polynomial with respect to 
E. The roots of this polynomial are spectral edges, 

2£+l 

W\E)=-^Yi{E-e,). (6) 
i=i 

The Bloch eigenfunction '5-1- and wavevector k are 
holomorphic functions with respect to complex vari- 
able E and they are defined on the Riemann surface 
^-W'^{E)^ This surface will be denoted by F. 

The A(a;, E) is a product of the two Bloch eigenfunc- 
tions of the Schrodinger equation ([3]). Since the eigen- 
function '^±(x,E) '-^ exjp(±i\^x) at £' — > 00, it follows 
that A{x,E) with respect to the variable £' is a polyno- 
mial of the ^-th order, 

e 

A{x,E) ^^+^_ ^CYl ilj {ix + uj)~E) 
i=i 

e 

= J2c,{E){p{tx + uj))'-^. (7) 

The roots ^n{ix -\- uj) are disposed inside gaps or on 
their edges and there is only one root for every gap. From 
the A(a;, E) definition, it follows that functions ^n{ix-\-uj) 
are periodic with a period 2|a;'|. When the variable x is 
changed, the function ^n{ix -\- uj) takes on values in the 
ri-th gap and extreme values of this function, which are 
defined by equation 9^7„(ia:: -I- w), coincide with the n-th 
gap edges. 
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In every band the A(a;, E) maintains the sign. When 
it passes from one band to another the sign changes. In 
the first band the A{x, E) > 0, Va;. It is easy to verify 
that A{x, E) satisfies the third order equation 

d^A4:[e{£ + f )p {ix + Lu)+ E]d.^A 

+ 2£{e+l)p^{ix + Lu)A = 0. (8) 

Functions ^j{ix + lu) or coefficients Cj{E) can be found 
if we insert the A{x,E) into the equation 

If we know function A(x,E), we can find the band 
edges by equation^ 

W'^iE) = -4[E + £{£ + l)p{ix + Lu)]A'^{x, E) 

- 2A{x, E)dlA{x, E) + id,A{x, E)f. (9) 

The wavevector k{E) is defined by expression 



k{E) 



where 



a;o + |2u;'| 

W{E) f dx _ W{E) I I 

4a;' 



A{x,E) 2i \A{x,E)/' 

(10) 



xo+2\uj'\ 



(/> = 



2U 



f{x)dx 



is the averaging with respect to the variable x. 

When E takes on values in a gap, the function A(2;, E) 
alternates its sign and the integral pil)) exists only in 
sense of the principal value. 

The function A{x,E) is periodic with respect to the 
variable x and has a period \2ui'\. It means that 



W{E) 
4w' 



dx 



A{x,E) 2uj 



- [F{x,E)+ik{E)x]+C, 



where F{x^ E) is a certain periodic function with respect 
to the variable x with a period \2u}'\. This function and 
the function under the integral have poles in the same 
place. 

The expression (llOp can be transformed into a more 
convenient form. We will use the following designation 



x{x,E) 



A(x,E) 
W{E) ' 



Then, if we pass to the new function xi^^ ^) the ex- 
pression ^ and differentiate it with respect to the vari- 



able E, we obtain 

dx" 



dx 



-2x-mi+^)pi^^ + ^) + E]-gj^ 

3x dx' 
X ' dE x dE' 



(11) 



Here and further, we denote a derivative with respect to 
X by prime. The function 

dx' X' dx 
dE x' dE 
is periodic with respect to x, and 



dx 

dE x' dE 



(12) 



Using formulas ^ and (fTT|) we obtain the following 
equality. 



Using formulas pOj) . p2|) and pS]) we derive the expres- 
sion for wavevector 



k[E) = ko + 



{A{x, z)) dz 



(14) 



The expression (fT4| defines a dependence of the elec- 
tron energy E on the wavevector k. This function is 
monotonous, but piecewise analytical, it corresponds to 
the extended band scheme. In the contracted band 
scheme the energy is periodic analytical function of 
wavevector in dual lattice with a period of this lattice 
(in our case this period is equal to 2'k/2\lj'\). 

Let us introduce £ parameters {tj} by means of the 
expression 



A{x,h,-,h) = ^+^-=C\{{p{ix + uo)-p{tj)). (15) 

Then, in terms of these parameters, we can present the 
functions \E'±(a;,i?) and k{E) in another form^^ 



\E'+(x, k) = u{x, k) exp{ikx) = 



a{tj — ix — Lu) 



a{ix + Lu)a{tj) 



exp j ix'^^ I exp(ifca;). 



(16) 
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'^-{x, k) = u{x, k) exp{—ikx) = JJ^ 



IX - 



a{ix + uj)(7{tj) 



exp j — ia;^^ j exp{~ikx), 



(17) 



/c(<i,...,i,) = ^(Cfe)-;7t, 



(18) 



r 



Inserting the expression lfT5|) into ([8]), we derive a sys- 
tem of the algebraic equationi^ to determine the quanti- 
ties p{tj): 

UkESk-i + Pk92Sk-2 + lkgzSk-3 - 5kSk = 0, k = l,£ 

(19) 

where 

ak=A{-lf-'ik~e~l), 

(3k = \{-lf-\2{£ - (fc - 4))3 - 15(^ - (fc - 4))2 
-f 37(^-fc + 4)-30), 

7fc = (-1)'= (£ - fc + 3) (£ - A: + 2) (^ - fc + 1) , 
4 = 2fc (-I)'' (4£2 + 2£ (2 - 3fc) + 2fc2 - 3fc + l) , 

and S*/!; is the symmetric function 



Sk 



^ p(iri)--p(irj- 



(20) 



'•l<r2<...<rfc 



Using the first equation in the system (|19|) . we can ex- 
press E through the quantities p{tj) and derive the fol- 
lowing expression 



(2£-l)^p(t,) 



(21) 



Below we describe in detail the spectrum of the 
Schrodinger equation with the one-gap and two-gap 
Lame potentials. 



A. The one-gap Lame potential {I = 1) 



The band edges are the following 

£l = 63, £2 = 62, £3 = 61. 



(22) 



The expression A(a;, E) for the one-gap potential is the 
following 

A(a;, E) = p{ix + to) - E = p{ix + oj) - p{t). (23) 

The spectrum of the Schrodinger equation ([3]) with the 
one-gap Lame potential in parametric form is the follow- 
ing 



E = p{t). 



(24) 



The expression for dependence k on E can be written in 
the form of elliptic integral 



k — ko + 



z + (ry'/c.') 



Bo 



2^(z - ei) [z - £2) (z - £3) 



zC?Z. 



(25) 



Lr Fig. [1] the complex plane and the fundamental do- 
main of the Weierstrass function ACHE are represented. 
When E takes on values in the first band the parameter t 
takes on values in the line AC. The point B corresponds 
to the minimum energy value, the points A and C cor- 
respond to the maximum energy value. When E takes 
on values in the second band the parameter t takes on 
values in the line DE. The points D and E correspond 
to the minimum energy value, the point corresponds to 
the maximum energy value (E — 00). 
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FIG. 1: The Weierstrass function fundamental domain on 
complex plane. 

Let us examine limit cases. 

In the limit case of free electrons the gap in the spec- 
trum vanishes. Let us put £2 = £3 = — iv' /'^') then 
£1=2 {r]'/u}') and 



k = 



dz 



2V(^-£i) 



= V(S-£i) 



or £^ = £1 + k'^. 

In the limit case of a very wide gap, we pass to another 
variable z = (1/2) • (£1 + £2) + (1/2) • (£2 — £1) in the 
integral ([25]) . Using expression 61 4- 62 + 63 = and 
introducing the designation A = {2E — £1 — £2)/(£2 — 
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£i), A £ we derive the following 



k = 



^ jsi + £2) + rf /lu' + ^ {62 ~ £i) if 



dip. 



Expanding the expression under the integral to the order 
of £2 — £1 and restricting ourselves to the zero approxi- 
mation, we obtain the spectrum in the limit case of tight 
binding electrons. 



where 



E ^ ^{£1 + £2) ~ ^ {£2 - £1) cos {ka) , 



2^3 



U^i+e2) + ^ 

It is possible to determine effective masses at band 
edges. According to the expression (pSj) . we have 



1 



y(e2-ei)(£3-£i) 2Vz-£i 
+ o((z-£i)i/2))dz, 



that leads to following expansions, 

^ + 0(fc3); _2(£2-ei)(£3-ei) 
ml 



E = £1 



ml 

fc2 1 



(ei +ri'/uj')'^ 



(26) 



E = £2 + — + o{k^y, — 



m 



2(£2 - gi)(£3 - £2) 

{£2+rj'/L0'f 



(27) 

Here we denote by m\ and m2 respectively the effective 
masses at the bottom edge and at the top edge of the 
first band. The effective mass m\ is positive and 
is negative. One can see that the both effective masses 
depend on all the band edges. 



B. The two-gap Lame potential (I = 2) 

For the two-gap Lame potential the band edges are the 
following. 



£1 



! = -3ei, £3 = -3e2, 

£4 = -363, £5 = yjig^. (28) 
The expression A(.x, E) in this case is presented in such 



a way. 



k{x, E) = l%p^{ix + uj)- 6Ep{ix + uj) + 2E^ 



-.92 



= 18[p{ix + cj) - p{ti)] [p{ix + cj) - p{t2)] 
'3 
2 



= 2 



(^p{ix + uj)- \Jg2- 3p2(ia;-f w)) 
{p{ix + cj) + \J g2 — 3p2(ia; -t- tj)j - i? 



(29) 



The spectrum of the Schrodinger equation ^ with the 
two-gap Lame potential in parametric form is the follow- 



ing 



= (C(tl)- ^^lj + (C (t2)-^t2), 



(30) 



The expression for dependence k on E can be written in 
the form of hyperelliptic integral 



3^2 



y/{z^ - 352) (8z 3 - 1832^ + 5433) 



dz, (31) 



When E takes on values in the first band the param- 
eters ti and t2 take on values in the Hue AC (Fig. [T]). 
When E takes on values in the second band the parame- 
ter ti takes on values in the line AC and the parameter t2 
takes on values in the line DE. In both cases p{ti) and 
p{t2) are real. The situation changes when E takes on 
values in the third band. In this case p{ti) and p(i2) are 
complex conjugate (fc and E remain real), ti and t2 take 
on values in the curves which are depicted inside fun- 
damental domain ACHE (while constructing, we used 
UJ = 1, Lu' — i). 

With the help of the Hcrmite reduction, we can express 
the integral ((3T|) through the standard elliptic integrals 
of the first and the second kinds^i^ and represent the 
spectrum in other parametric form. Let us put 



2z' 



y 



3(z2 
Using equalities 

E 

zdz 



ay 



3ff2, h = -54^3. 



^{z^-a) {Sz^-Qaz-h) 2V3J y/y3~3ay + b' 

Eq Yq 



dy 



(2z2 _ a)dz 



^{z'^-a) {%z^-Qaz-h) 



]z^ — 6az — b 



ydy 



2V3J ^y3^3ay + b' 

and putting y = 6p(h), we can write the integral pip in 
terms of elliptic functions 



E'\+27g3 
9B^-27g2 



.(32) 

In this form the spectrum of the two-gap potential is 
provided in the paper 
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By the same way as in the case of one-gap potential we 
derive expressions for effective masses in the first band, 



I 



2(e2 - £i)(£3 - ei)(e4 - £i)(£5 - £i) 



((£i)2+3(77Vw')£i 



552 



E = e2 



I 



2(e2 - ei)(£3 - £2)(e4 - £2)(£5 



((£2)2 + 3(r?'/a;')£2 



(33) 



(34) 



252 



where m\ and are effective masses at the bottom and 
at the top band edges respectively. One can see that the 
effective masses depend on all the band edges, the same 
as in an one-gap case. In Table HI the effective mass 
values are provided for different values of the parameter 
|r|, we put a period of the potential T = \2uj'\ = 2. 

TABLE I: Inverse effective masses for the first band in case 
of the two-gap Lame potential. 



1/ral 



0.55 
0.74 
1 

1.23 
1.52 
2.18 



1.99 
1.89 
1.25 
0.57 
0.15 
0.005 



-47.34 
-10.03 
-2.13 
-0.67 
-0.16 
-0.005 



In Table HI one can see that with the increase of the 
value of |t|, the values of |1/»ti5;| and |f/?7i2| decrease. 
When the band degenerates into an energy level that 
corresponds to the limit case of tight binding electrons, 
TO* 00, 7712 ^ —00. We obtain so called "heavy" elec- 
trons. In the other limit case of free electrons, 777^ tends 
to the mass of free electron, i.e. m\ 1/2. And since 
the gap width vanishes, TO2 — > —0 in this case. 

In Fig[2l we provide the graphs E — E{k) and U{x) for 
the two-gap potential corresponding to different values of 
the parameter |t| = \lo' /lo\. As |t| increases, the width 
of the gaps increase. There are the limit cases of very 
narrow gaps and very wide gaps in this figure. With 
increase of |t|, the second band also degenerates into an 
energy level. 
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FIG. 2: Dispersion relation E = E{k) and potential U(x) in 
case of the two-gap Lame potential for different values of the 
parameter |t| = \uj' /uo\. 



III. THE SEPARABLE FINITE-GAP LAME 
POTENTIAL 

The one-electron separable Hamiltonian has the form 

+ U,{x,)\. (35) 



dx 
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In general case n € N, but we consider below only the 
three dimensional case, n = 3. We will use the finite-gap 
Lame potentials for the potentials Uj{xj). In general 
case of the orthorhombic lattice every one-dimensional 
potential Uj(xj) is characterized by its own parameters 

H = -dl~ dl - dl - 4(4 + i)p{ix + c^,) 

- -h \)p{iy + UJy) - 4(4 + ^)P{IZ + UJy). (36) 

In Fig. [3l the three dimensional one-gap potential for 
z = is represented with the following half-periods: lo^ = 

UJy ^ I, U}'^ = Uj'y = i. 




FIG. 3: The three dimensional one-gap Lame potential for 
fixed value of the variable z, we put z — 0, ujx = ujy = 1, u)'^ = 

An eigenfunction of the Hamiltonian (|36p is a product 
of one-dimensional eigenfunctions of every Hj , 

3 

*i,k(x) = n^^-'=.(^^)' (37) 

where I (4, 4, 4), Ij £ {1,..,4} is a band number. 
The expression for energy consists of three terms, 

Ei(k) ^ Ei,{h) + Ei.ih) + Ei.ih). (38) 

The Fermi surface is described by the following equation 

Ef £;i(k + b), 

were Ep is the Fermi energy. Here the wavevector k 
takes on values in the Brillouin zone, every solution of 
this equation for a prescribed value of dual lattice vector 
b defines one of sheets of the Fermi surface. The set of 
dual lattice vectors corresponds to their set in Fourier 
series expansion of a separable potential. If the real elec- 
tron potential differs a bit from the separable one, then 
the corresponding corrections can be taken into account 
using the perturbation theory. 

Let us examine the limit case of very narrow gaps. 
Since the potential periods Tj = Sjcj^l are fixed, let 
ujj — > oo. Then, potentials Uj{x) tend to constants, the 



dependence of energy on wavevector becomes quadratic 
and the Fermi surface sheets are similar to pieces of 
spherical surface. Thus, our model is the considerable 
generalization of the Harrison method that was success- 
fully used to construct the Fermi surfaces of metals jiS, 

Let us examine the other limit case of very wide gaps. 
Since we fixed the potential periods Tj = 2|a;^|, let 

LUj — > 0. Then, — eg"''' —i- 0. It is obvious that 
in this case all the bands degenerate into levels. The 
finite-gap spectrum of tight binding electrons is similar to 
the spectrum that is obtained in the well-known LCAO 
method. In this case function 'i'l,,^ (x) ; 4 = (0,0,0) 
corresponds to s-state. The three functions ^'/j (x) , 
^i,,o{x), "Si^oix); 4 = (1,0,0), 4 = (0,1,0)', 4 = 
(0, 0, 1) correspond to p-state as they are proportional 
to X, y, z respectively when the space variables tend to 
zero. The six functions "^j^^q{x), '^i^^o{x), '^jg olx), 
^i,, o{x), o{x), o{x);'l4 = (1, 1, b), 4 = (1, 0, 1), 
4 = (0,1,1), 4 = (2,0,0), 4 = (0,2,0), 4 = (0,0,2) 
split into five functions that correspond to d-state and 
one function corresponds to the s-state since these func- 
tions are proportional to a;?/, xz, yz,x'^ — a, y^ — b, z^ — c 
respectively when the space variables tend to zeroi^ 

The Fermi surfaces for cubic lattice are represented in 
Fig. [3][5l They are constructed using the separable one- 
gap Lame potential for different values of the parameter 
|t|. In Fig. [31 the Fermi surfaces in the first band are 
represented at Ep — 2.43. The Fermi surface sheets don't 
appear in the second band at this value of Ep. In Fig.|4l 
the open Fermi surfaces are represented at Ep = 3.8 in 
the first band. In Fig. O the closed Fermi surfaces are 
represented at the same value of Ep in the second band. 

IV. THE LIFSHITS SINGULARITIES OF THE 
ELECTRON THERMODYNAMIC POTENTIAL 
AND THE ONE-GAP LAME POTENTIAL 

It is known that the density of states in a crystal 

dN{E) _ 2V f dS 
^ ' dE (27r)3 J |V£;(k)|' 

as a function of the variable E has singularities which are 
named the Van Hove onesi^ They occur when the elec- 
tron group velocity vanishes, Vi?(kc) = 0. Since i?(k) 
is a smooth function of the wavevector k, it can be ex- 
panded into the Taylor series about a critical wavevector 
kc. Restricting ourselves to quadratic terms, we obtain 

where ^ — k — kc. Depending on signs of the coeffi- 
cients ai, a2, as the four types of critical points are dis- 
tinguished, they are a minimum value, a maximum value 
and two types of saddle points. In a three dimensional 
crystal the function i^iE) at the critical energy Ec is con- 
tinuous but its derivative at this point has infinite dis- 
continuity. 
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\T\=1.2i 



FIG. 4: The Fermi surface in the first band for different values 
of the parameter |t| = /uj\, the Fermi energy is Ef ~ 2.43. 




I T 1= 0.73 



|r|=l 



FIG. 5: The Fermi surface in the first band for different values 
of the parameter |r| — /uj\, the Fermi energy is Ef = 3.8. 



The function i?(k) is periodic with respect to every 
component of the wavevector k — (k^, ky, k^), it means 
the energy Eikx^ky^kz) is defined on three dimensional 
torus. This function has a denumerable set of branches 




FIG. 6: The Fermi surface in the second band for different 
values of the parameter |r| = \uj' the Fermi energy is 
Ef = 3.8. 



and every one of them corresponds to a certain energy 
band. The function E{kx,ky,kz) is difFerentiable and 
regular so, according to elementary facts of topology, for 
each of the branches the number of critical points re- 
lates to the Betti numbers of the manifold, where the 
function E is being defined. Thus, every branch of the 
function E{kx, ky, kz) has one maximum value, one min- 
imum value and three saddle points of every type. But 
in a real crystal the bands may overlap and as a result 
the singularities of the function v{E), that are related to 
different branches, may compensate each other. 

In 1952 Van Hove examined singularities of the elas- 
tic frequency distribution of a crystal (that are similar to 
singularities of the electron density of states) and showed 
that for a three dimensional crystal the frequency distri- 
bution function has at least one saddle point of every 
types and its derivative takes on the value — oo at the 
upper end of the spectrum. '^^ 

In 1960 I.M. Lifshits connected singularities of the den- 
sity of electron states with a variation of the Fermi sur- 
face topology,^*' He showed that at zero temperature the 
electron phase transitions occur at the special energies, 
they were later named the Lifshits electron phase transi- 
tions of 2i kind. Singularities of the conduction electron 
spectrum at the point of such a transition lead to anoma- 
lies of thermodynamic and kinetic quantities. When tem- 
perature increases, this transition smoothes out. 

At present, the Van Hove singularities are considered 
as one of the reasons for the origin of high-temperature 
superconductivityi^i 

Let us consider the Schrodinger equation with the one- 
gap separable Lame potential that corresponds to the 
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FIG. 7: The Fermi surface in the first band in the form of a 
goffered cylinder, we put uj = 1, uj' = 0.65i; u = 0.3, u' = 
0.5i. 



tetragonal Bravais lattice, 

-A* + W{x, y, z)* = 
W{x, y,z) = — 2p{ix + Cj\Cj, a)' ) 
—2p{iy + (Dp, a)' ) — 2p(iz + wjo;, cj' ) — 63 + 4ei. 

In order to describe the Lifshits electron transformation, 
it is enough to examine the case when the Fermi surface is 
open only at two opposite sides of the Brillouin zone. Let 
us choose half periods and the Fermi energy so that the 
Fermi surface has the form of a goffered cylinder which is 
prolate along the direction , as in Fig. [71 In this case 
we have to put a) > | > \lu' \ (while constructing, we 

used UJ = I, uj' = 0.65i; uj — 0.3, uj' — 0.5i). Then, in 



order to simplify computations, we can expand potential 
with respect to the variables x, y and restrict ourselves 
to zero approximation. 

The spectrum of the Schrodinger operator with the 
potential 



U{x,y,z) 

looks as follows: 

E[kx, ky, hz) 



-63 



2p{iz + uj\uj, uj' ) 



Cihz) ~hhz 



(39) 



(40) 



The function E{kx, ky, kz) has two critical points in the 
first band that are a minimum value and a saddle point 
and one critical point in the second band that is a mini- 
mum value. We will consider below only the first band. 

One can see that -E|k=o — 0. Let us put K'^ = kl + ky. 
For the first band hz — t + uj' , t G [—uj; uj]. Since, the 
function E(k) is even, it is enough to consider the case 
when t e [0; uj]. The expression for the number of states 
is the following 



N{E) = 
2V 



2V 



(27r)3 



2V 
(2^ 



nK^dkz 



(27r)3 
2V 



TT{E + e3~ p{hz))dkz 



(2 



^ j 7r{E + es-p (hz)) (^-p {hz) - ^ j dhz 



(41) 



where V is dimensionless volume. 

It easy to show that at the energy Ec = —e^ + 62 the 
Fermi surface opens on the boundary of the Brillouin 
zone which is orthogonal to the axis kz. In this case we 
can express the number of states (j41l) in such a form. 



to 

N{E)^^ J {E + e3-p{t + uj'))(^-p{t + u;')-^^dt 



V 

2^ 



^-'±iE + e,)]to + {ato+uj')-7^')UE + es)-^, 



+ -p' {to + uj') 



(42) 



r 



Here we have introduced the parameter to as follows. B. Under condition E > Ec we define the parameter 

A. Under condition < E < E^ we define the param- to by the equality tp = ^■ 
eter to by the equality 

E + e3^p{to + uj') . 

According to this equality, the function to ~ to{E) is 
differentiable and has a derivative 

dto _ 1 

dE ~ p' {to + uj') ■ 



Differentiating this expression for the number of states 
in E and using the definition of parameter to we obtain 
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the electron density of states, 



27r2 



2^(C(io + -')-^(to+c.') 



^ 1 E>E, 



2t! 2|u;'| ■ 



(43) 



Even though this expression was derived for the one-gap 
potential, it is correct in a case of any £— gap potential. 

Let us expand the function v{Ei) under the condition 
E < Ec \vl the ^Te powers. We derive the following 
formula, 



v{E) = 



V E^'^ , ,,,, 6e2-f 



27r2Ci/2 



2(e3 + ^)^ 



(44) 



It is obvious that the first term in the expansion 
corresponds to the function E{kx,ky,kz) — k'^ + ky + 

Let us calculate the first derivative of the density of 
states. We should distinguish two cases. 

A. Under condition < E < Ec we obtain the expres- 
sion 



diy{E) 



dE 27r2 p' (io -I- uj' ) 

V & + e, + E 



47r2 £;i/2(ei-e3-£;)i/2 (^^ _ £;)i/2 



(45) 



B. Under condition E > E^ wc obtain the expression 
dv{E) 



dE 



0. 



(46) 



The singularity at the point E — ei — corresponds to 
the Fermi surface appearance in the second band. The 
singularity at the point E — Q corresponds to the min- 
imum value of E in the first band (the Fermi surface 
appears in the center of the Brillouin zone). The sin- 
gularity at the point E — Ec corresponds to the saddle 
point (the Fermi surface opens on the boundary of the 
Brillouin zone). In all three cases the singularities have 
positive signs. We will consider further only the case 
-E = i^c, it is obvious that the other two cases can be 
considered similarly. 



Let us introduce the following notation: 

n' 



a 
2 



V 



62 



(62-63)1/2(61-62)1/2 

l^ + e2| 



V 

4^' (^-63)1/2(61-62)1/2' 
TT rjrj' 



C 



— •( 


'92 


27r2 \ 


L 

,12 


V 


1 


2tt 2 




V 


\el 


"2^ 


2 




\uj' 12 



UJ 



TT 



2\uj' 

1 



, 92 
62 1 



7^? 



40 



52?7 



33 



20 

In order to compute the thermodynamics poten- 
tial il(T, y,/i) we expand the following expressions 
dv{E)/dE, v{E), N{E) in the {Ec - Ef''^ powers. We 
derive the following expressions: 

( a 1 

dv{E) 



dE 



uiE) = 



N{E) = 



2 (_E^-_E)i/2 

+0{{Ec-Ey/^), E<Ec, (47) 
0, E> Ec, 




3/2^ 




E) 



BE + la{Ec - Efl"^ 
+6 {{Ec - Efl^) 
BE, 



E<Ec, (48) 
E>Ec, 



E<Ec, (49) 
E>Ec. 



The coefficient a is expressed in terms of the effective 
mass 7712 along the k^ direction at the top edge of the 
first band, 

V 



\/27r' 



rrin 



(50) 



Hence the coefficient a depends on all the band edges and 
a e [0;oo). When the gap width decreases on boundary 
of the Brillouin zone, which is orthogonal to the kz axis, 
the value of the coefficient a decreases also. In the limit 
case of extremely wide gap the parameter a —^ 00. In this 
case the problem from a three-dimensional one becomes 
two-dimensional and we should use other speculations. 
In the other limit case when the gap width vanishes the 
parameter a —^ 0. It is obvious: when the gap width 
equals zero the value Ec is not a critical point. 

Thermodynamic potential at a normal metal is defined 
by the following expression 



n{T,v,fi)^^ 



N{e)de 



1 + exp (^) ' 

u 

n{0, V,fi) = no = - J N{e)de 



(51) 
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Using the expression (|42|) we can calculate the integral 
([5T|) and find the analytical expression for S1(0, V, /i), but 
this expression will be quite a complicated one and we 
do not adduce it. 



Using the formula (|49|) we present the electron thermo- 
dynamic potential Qq in a neighborhood of the critical 
value of energy in such a form, 



J 



^ r C + AE, + ^BEl - Afi - ^B^i'' + ±a[E, - y.fl^ + O {{E, - ^.y'^) , /i < E,, 
° \C + AEc + \beI- Aii-\B^i'^, pi>Ec. 



(52) 



The variation SN{E) of number of states looks as fol- 
lows: 

SN(E)^l -3°^^^^-^^'^'' ^<^- 
^ ' 1 0, E>E^. 



5N{e)de 
1 + exp {^) 



Then the variation of the electron thermodynamic poten- We provide below the asymptotic expansion of the ex- 



tial can be written down in such a way 



pression 6fl(T, V, /i) under condition T ^ 0, 



sn 



-^aET^ + ^a{E, - ^ifr2 + ^^^T^E, - /i) 



5/2 



U/2 



|gc-Ml 



0{T^ 
On exp 



T 



(53) 



A similar expression was presented in the paper 

In Fig. [51 the graph d'^Q,{T, V, /i) on the temperature 
T and the chemical potential fi is represented (while con- 
structing, we used V — 20, uj ~ 0.3, uj' ~ 0.5i, accord- 
ingly Ec = 2.5). One can see that when /i Ec, T = 
the third derivative from the thermodynamics potential 
has an infinite discontinuity and 



lim dln{0, V, /i) = -oo 



One can see also that the singularity at the point 
fj, = Ec smoothes out with the temperature increase. 
In order to estimate an appropriate temperature inter- 
val, we should pass from dimensionless units to Kelvin 
degrees. In order to do this we have to multiply the di- 
mensionless temperature by the following coefficient. 



Eo ^ 

kh 2ma?'k[, 



2m(c/2|w' \ykb 



(54) 



where fc^ is the Boltzmann constant and c is a period 
of crystal potential. If the potential period equals iA 
and effective electron mass is put lOrrie, where me - 
mass of free electron, we find out that the dimension- 
less temperature value 0.1 corresponds to the tempera- 
ture value in standard units ^ A%'^K. Since the width of 
the area A{Ec — /i), which contributes to the singular- 
ity of d'^fl{T,V, II), is proportional T, then in a general 




FIG. 8: Third derivative of thermodynamics potential, we 
put V = 20, u ^ 0.3, uj' = 0.5i and = 2.5. Thermody- 
namics potential, temperature and chemical potential take on 
dimensionless values. df^QlT, V, fx) has an infinite discontinu- 
ity when n — > Ec, T = 0. 



case the singularity smoothes out slowly enough with the 
temperature increase. 
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V. MAGNETIZATION SINGULARITIES OF where /3 = eh/2mc is the Bohr magneton, H is an exter- 

THE PAULI PARAMAGNETICS nal magnetic field, T is a temperature. According to this 

expression the magnetization at T = can be presented 
A magnetization of an electron gas in magnetic field is iii the following way, 
described by the following expression. 




M(T,n,H) = - / I r r B 



(56) 

, , Let us expand the expression ([5^)1 in the (E^ — fi±l3HV'^ 
-1 I / s-ii+i3h \ I ^' ^ ' powers. Using the formula ([49]) we come to the following 



' ^{2PHB + la{E, ^n- /3i7)3/2 „ 2^(^^ _ ^ ^ ^^)3/2 

+0((S,-/i±/3ff)5/2)), ti<E,-f3H, 

§{2f3HB - la{Ec - + PH'f/'^ + 0{{E^ - + PHf^)), E^ ~ /3H < < E^ + f3H, 

^ P^HB, fi> Ec + PH. 



Since d^u{^± PH) has infinite discontinuity when /i± 
PH = Ec then the second derivative of the appropriate 
magnetic moment with respect to the chemical potential, 



d^M{0, p, H), has singularity also when /i ± PH — Ec- 
The sign of singularity is defined by the sign ahead of 
N(fi ± PH) in the formula for magnetization. 



= \ - 4(g.-;+%i/2 + 0{{Ec - + PH)'/^), Ec-pH <fi<Ec + PH, 

0, " H> Ec + PH. 



The expression (|50p defines the relation of the coefficient 
a to the effective mass As mentioned above, with the 
decrease of gap width on boundary of the Brillouin zone 
a value of the coefficient a decreases. In the limit case 
when the gap width equals zero the coefficient a = 0. In 
this case the value Ec is not critical point and there is no 
singularities that is related to this point. 

In Fig. [21 the graph M(T, fi, H) on the temperature 
and the chemical potential is represented (we put PH = 
0.15£;c)- 

In Fig. [ini the graph d^M{T, fi, H) on the tempera- 
ture and the chemical potential is represented (while con- 
structing, we used PH = 0.15i?c, w = 0.3, uj' = 0.5i, ac- 
cordingly Ec = 2.5). One can see that when ^ Ec^PH 
at T = the second derivative from the magnetic mo- 
ment has infinite discontinuities and 



lim 



. dlMo 



oo, lim d,,Mo — const, 



lim 



. dlMo 



lim 



dl^Mo — const. 



As in the previous case, the singularities at the points 
Ec^PH smooth out slowly enough when the temperature 
increases. 



VI. CONCLUSIONS 

In this paper, we have examined the one-gap and two- 
gap separable Lame potentials in detail. We have con- 
structed the dispersion relation E(k) for one-dimensional 
case and the Fermi surfaces in the first and second bands 
for three-dimensional case. The pictures illustrate a pas- 
sage from the limit case of free electrons to the limit case 
of tight binding electrons. We have provided also ana- 
lytical expressions for effective masses of electrons in a 
metal. These expressions depend on all the band edges. 

These results have been used to study the singulari- 
ties of the electron part of the thermodynamic charac- 



13 




teristics in metals. We have derived explicit analytical 
expressions for the density of states in a metal and its 
derivative. Then we have examined the thermodynamic 
potential and magnetic moment of metal in an external 
magnetic field and extracted the singular part. We have 
also obtained the relation between the parameter of sin- 
gularity and corresponding effective mass. All the expres- 
sions we derived contain the parameters of our potential 
only. Therefore, if we fix potential, we can calculate all 
the coefficients of these expressions. 



FIG. 9: Magnetization of an electron gas in magnetic field, 
we put I3H = 0.15Ec, to = 0.3, to' = 0.5i and Ec = 2.5. 
Temperature and chemical potential take on dimensionless 
values, magnetization is measured in Bohr magnetons. 




FIG. 10: Second derivative of electron gas magnetization in 



magnetic field, we put PH — 0.15i5c 



0.3, u = 0.5i 



and Ec = 2.5. Temperature and chemical potential take 
on dimensionless values, magnetization is measured in Bohr 
magnetons. d^M{T, /i, H) has infinite discontinuities when 



Ec±l3H,T = 0. 



It is necessary to point out that instead of electrons, 
the arbitrary quasi-particles can be considered within a 
framework of this approach. 



We would like also to say a few words about generaliza- 
tion of the finite-gap potentials to lattices of other spacial 
symmetry. In the paper— the authors have tried to de- 
velop such a generalization for HCP lattice by means of 
the perturbation theory and calculated the Fermi surface 
of Beryllium. The elliptic Calogero-Moser potentials are 
also of significant interest concerning their use in solid 
state physics. These potentials are multidimensional. 
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